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1 INTRODUCTION 



Hadroproduction of the particles with open charm, among other reactions 
(e + e~ — > cc etc) is of interesting possibility to study a hadronization process 
of the heavy quarks, such as c — > D, A c , ... (]]]. 

It is well known that the production of charmed particles in e + e~-anni- 
hilation, 

e + e- -4 DI, (1) 

may be separated into a process of cc-pair production and a process of inde- 
pendent fragmentation of each c (c) quark into .D-meson: 

c(c) -> D(D)X. (2) 

In so doing the fragmentation process is described by means of fragmenta- 
tion function D(z), where z = Pd/Pc {pd and p c are the momenta of .D-meson 
and c-quark, respectively). 

It worth to note that a scale-invariant description becomes available at 
sufficiently high energies, i.e. in the limit of ^/s ^> mo- For the low energy 
case ( \fs > tub) the role of non-scaling (power) corrections is very essential. 
As a result at that energy region one can not represent the cross section in 
a simple factorized form: 

^^a c - c ®D{z). (3) 

The question arises of what extent a mechanism of fragmentation of heavy 
quarks, produced in e + e~ annihilation, can be applied for the case of hadronic 
interactions ? 

Note, that in hadronic interactions the presence of light partons from the 
initial hadrons will affect the hadronization scenario. One should distinguish 
two kinematical regions, namely, the region of p-\ < po and the region of 
Pt > Po- Like in e + e~ annihilation, the fragmentation process above po has 
a fragmentation character and the cross section has the form as follows: 

dan r da, 



®D{z)dz. (4) 

dpj J dk-j 

At the region of pj < po the presence of light partons, produced simultane- 
ously with cc pair, may radically altered the factorization in the form of (13). 



2 



In our previous publication [0 we made an effort to take into account 
this fact by consideration of a some contribution into hadronization as a 
process of recombination with valence quarks. The rest part is described by 
a fragmentation mechanism. In the framework of such a consideration one 
can properly describe the asymmetry in the production of leading and non- 
leading charmed hadrons 0. Note, that our description of total spectrum 
is not well adequate. It explains because the excess of our calculated distri- 
bution at low x. This fact is caused by the application of the fragmentation 
model for the whole kinematical region. 

In the present article we make an effort to improve the description of c- 
quark hadronization at low pt region. We try do not use the fragmentation 
model as a basic mechanism. We apply this model only at boundaries of 
the phase space, where c-quark momenta achieve their maximum values and 
where the hadronic accompaniment of c-quarks is very small. 

We note also, that in the present article we do not consider the question 
related to an absolute magnitude of the cross section of charmed particle 
production. The matter is that such a problem is determined by taking into 
account the higher order corrections of the perturbative QCD as well as the 
choice of the strong coupling constant of a s (/i 2 ). Indeed, the analysis of 
0(cel) corrections to cc pair production cross section (see, for example, 0) 
exhibits that their inclusion practically do not alter the form of the inclu- 
sive distributions of c-quarks. Therefore, in what follows we shall restrict 
our consideration to study the form of differential spectra of the charmed 
particles. 

The article is organized as follows. We consider the modification of the 
fragmentation mechanism in the Section 2. We compare in Section 3 the re- 
sults of our calculations with the experimental data on .D-meson production 
in 7r-beam. Section 4 presents the model predictions for the K beam. The 
main results are summarized in Conclusion. 
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2 c-QUARK HADRONIZATION IN HAD- 
RONIC COLLISIONS 



In the parton model, the cross section for the production of heavy quarks Q 
in hadron-hadron collisions has the form 

a(hh 2 -> QQX) = £ M - J a(ij -> QQ)ft 1 (x 1 )dx 1 f} 2 (x 2 )dx 2 , (5) 

where summation is performed over all partons participating in Q quark pro- 
duction ij — > QQ, a is the cross section for the corresponding hard subpro- 
cess, fi(x) is the distribution of z-type partons in the hadron h, Xi( 2 ) is the 
fraction of the initial-hadron momentum carried away by the corresponding 
parton. 

As it mentioned in the Introduction, the presence of light partons from 
the initial hadrons leads to radically different scenarios of the hadronization 
of charmed c-quarks, produced in hadron-hadron collisions and in e + e~ an- 
nihilation. In particular, the interaction in final state with valence quarks 
from initial hadrons (recombination) makes possible to explain the leading 
particle effect in the charmed hadrons production [0, [!| (i.e. the difference 
in x-distributions of D and D mesons, as well as A c and A c baryons). 

Such an inclusion of interaction of charmed quarks with valence quarks 
from the initial hadrons is carried out by the introduction of the quark re- 
combination function of R(xy, z;x) 0. Below we present the basic concepts 
of the recombination mechanism. The detail description of this mechanism 
is given elsewhere [0, ^J. The recombination of qy and Q quarks into Mq 
meson is described by the function of R(xy, z; x): 

R(x v ,z;x) = p(£ v ,£q) - £v - 6q), (6) 
/£ t \ _ T{2- a v - a Q ) c (i- av ) c {i- aQ ) m 
PW,£q) ~ T(l- av )T(l-a Q )^ v > (7) 

where = xy/x and £q = z/x, while xy, z, and x are the fractions of the 
initial-hadron cm. momentum that are carried away by the valence quark, 
quark Q, and the meson Mq, respectively, ay and ccq are the intercepts of 
the leading Regge-trajectories for the qy and Q quarks, respectively. In our 
calculations we use 0. 

1 

a u = a d = -, a s ^0, a c « -2.2. (8) 
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With the aid of the function R(xy, z; x) describing the recombination of 
the quarks qy and Q into a meson, we represent the corresponding contribu- 
tion to the inclusive spectrum of Mq mesons as follows: 

x *~3 — = Rq I x vz* f °-, R(x v , z; x)^-—, (9) 
ax J dxydz Xy z 

where x* = lEj^fs and x = 2pi/yfs (here, E and pi are the energy and 
longitudinal momentum of the Mq meson in the c.m.s. of the initial hadrons); 
xy and z are the momentum fractions carried away by the valence quark and 
heavy antiquark, respectively, and xv z * d f v a dz is the double-differential cross 
section for the simultaneous production of the quarks qy and Q in a hadronic 
collision. 

The parameter, Rq, is the constant term of the model, that determines 
the relative contribution of recombination. In the present article the best 
description of the experimental data could be achieved at 

Ro « 0.8. (10) 

Note, that using of the recombination with the valence quarks is indis- 
pensable to an explanation of the leading particle effect. At the same time its 
contribution in the total inclusive cross section production of the charmed 
particles is sufficiently small (~ 10%). This mechanism dominates in the 
high x region. 

In the conventional fragmentation mechanism the inclusive cross section 
for the charmed hadrons (D mesons) production has the form as follows: 

T7i d 3 a F f dV(feife 2 -> cX) 3 

E H — — = / E c D(z)5{p H - zp c )dy c . (11) 

d 6 pu J d 6 p c 

In the low x region, which determines the basic contribution into the total 
cross section for charmed hadron production, the hadronization process of c 
quark has more complicated nature. 

Indeed, when evaluating the spectra of charmed particles produced in 
hadronic collisions one assumes that the fragmentation function D(z) is well 
known from other experiments (in particular, from e + e~ annihilation). One 
of the widely used parameterization is as follows : 

D(z)~[z(l-±--^-)}-*, (12) 
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where the parameter e ~ m 2 / m?Q is determined by the type of a hadron (for 
instance, one has e D o = 0.135 ± 0.010 and e D * = 0.078 ± 0.008 0). 

Another parameterization, proposed by us early |J, takes into account 
the Regge-asymptotic at z — > 0: 

D(z) ~ z- a Q{l-zy, (13) 

where 7 ~ 1, while aQ is the intercept of the leading Regge-trajectory for 
the Q quark (a c ~ —2.2, see (§)). 

Both of these parameterizations provide a reasonably fair description of 
the experimental data. The comparison of description of the reaction e + e _ — > 
-DX by means of the fragmentation functions in both form of ( 12]) and (|13|) 
is presented in || |£J] . 

As it mentioned in Introduction, the use of the fragmentation function is 
justified at asymptotically large values of the invariant mass of cc pair, namely 
at M c5 ^> 2m c . However, this condition is not obeyed in the hadronic pro- 
duction of the charmed particles. For such the principal contribution 
into inclusive charm production cross section results from c quarks with low 
values of the invariant mass of cc pair (M cS > 2m c ). These quarks dominate 
in the central region on Feynman variable of x. The pairs of c quarks with 
a large invariant mass, where one may apply the fragmentation formalism, 
give a noticeable contribution at high x and high pt, as well. 

These arguments are illustrated in Fig. 1. This picture presents the in- 
clusive x-distribution for c quark for all M c5 (the upper histogram) and for 
the c quarks with the invariant mass of M cE > M Q = 10 GeV (the lower 
histogram). As is seen from the figure just the charmed quarks with small 
invariant masses of cc pair give the dominant contribution in the charm 
production cross section in the central region, while the region of x — > 1 
corresponds to the contribution due to a large values of mass of M c5 . 

In this figure we also present the spectrum of charmed particles summed 
over all types D and D mesons (for the reaction of tt~N collisions at E n = 



250 GeV |10| .) This experimental spectrum should be compared with the 
distribution of c quarks. One may deduces from the Fig. 1 that the correspon- 
dence, like a duality relation, takes place. Namely, the spectrum of charmed 
hadrons, summed over all types of charmed mesons, is well described by the 
inclusive spectrum of c-quarks. 

Such a satisfactory description of summed spectra of D mesons by pure 



c-quark spectra was also pointed out early [p], |I0"|. However, it is evident that 
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in the framework of pure fragment at ional mechanism one should expect the 
identical spectra of both D and D mesons (as well as charmed baryons and 
antibaryons). As a result, one can not reproduce the leading particle effect. 
As it mentioned above, this effect can be described with the help of the 
recombination mechanism. Therefore, one could assume that the inclusive D 
meson production cross section is described by the sum of two mechanisms 
as follows: 

dap = da§ F (p D = p c ) + da% c 
dx dx dx 

where the first term corresponds to the "hard" (HF) fragmentation (in that 
mechanism a charmed quark does not lose its momentum in hadronization 
process), while the second term corresponds to recombination contribution 
(that mechanism takes into account the charmed c-quark interaction with 
valence quarks from initial hadrons). 

However, such a simple addition of the recombination contribution to the 
c quark spectrum (i.e. D meson) does not provide to reproduce the behavior 
of ^-dependence of the corresponding asymmetry A: 



^ -^(leading) — ^(non — leading) 
^(leading) + ^(non — leading) 



Indeed, the Fig. 2 presents the description of the asymmetry (|TjJ) by 
means of the equation (|T3). The different histograms in this figure correspond 



to the different values of the parameter R (see the equation (g) for the 
recombination mechanism). As is seen from this figure it is impossible to 
obtain the simultaneously proper description of the asymmetry A both at 
small x (0 < x < 0.4) and at high x (0.5 < x < 0.8) as well as with the 
description of the D mesons inclusive x-distributions. 

This is because the simple equating of D meson spectrum to the c quark 
spectrum (i.e. D(z) ~ 5(1 — z)) leads to exclusively "hard" spectra of D 
mesons at high x. Therefore, we get the conclusion that one needs to use 
a more softer (in compare with 5(1 — z)) fragmentation function for the 
description of charmed quark hadronization at high values of the Feynman 
variable of x. 

Bearing in mind the preceding, we consider the modification of the con- 
ventional fragmentation scenario of charmed quarks hadronization. It seem 
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likely that in the region of small invariant mass of cc pair the description of 
the hadronization in terms of fragmentation function (evaluating from the 
e + e~ annihilation) is not justified. Indeed, one has large amount of partons 
from initial hadrons in the central region of x. Therefore, the c quark in com- 
bination with one of such a parton can easily to produce a charmed hadron. 
Such a process occurs practically without any loss of c quark momentum (i.e. 
p D pa p c ). Therefore, in the small x region one should expect the coincidence 
of the spectra of D mesons and c quarks. Whereas at high x region one 
may use the conventional fragmentation mechanism (and a recombination 
mechanism, as well). 

Hence, we consider two regimes of the charmed quark fragmentation. 

a) Close to the threshold of c quarks production, i.e. at M c5 > 2m c , the 
momentum of the produced D meson should practically coincides with the 
momentum of the charmed parent-quark. 

b) For cc pair with the invariant mass M c5 greater that certain scale of 
Mo (where Mq 3> 2m c ) the c quark hadronization process may be described 
with the help of fragmentation function (for instance, in the form of (|T2]) 
or®). 

In terms of the fragmentation mechanism these two regimes may be repre- 
sented in a uniform way by introduction of the dependence of fragmentation 
function from an invariant mass of cc pair as follows: 

b mf, M \ j ~6{l-z) at M c - C « 2m c 

1 ' c5) \ D{z) from (12) or pD at M cc - > M 1 ' 

It worth to note, that assumed M cg -dependence of the c-quark fragmen- 
tation function has nothing to do with logarithmic violation of scaling in the 
fragmentation function. 

To reproduce in a uniform way the both two fragmentation regimes ([16|) 
we use the simplest expression for D(z, M c5 ) in the from of (0) as follows: 

D MF (z, M c5 ) ~ z - a{M ^\l - z), (17) 

with two additional conditions on a(M cS ): 

a(M c5 ) — > -oo at M c5 -> 2m c , , . 

a(M c5 ) -> a c at M cc - w M . 1 ) 
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Our parameterization for a(M c5 ) is presented in the Appendix 1. A fit to 
experimental data exhibits that the magnitude of the parameter Mo, 

M « lOGeV, 

provides a satisfactory description of experimental results. Such a magnitude 
does not contradict to the experiments in e + e~ annihilation, where already 
at the energy of y/s ~ 10 GeV one can describe the process in terms of 
fragmentation mechanism by using the dependence in the form of 

Thus, the summed differential cross section production of charmed hadron 
(D-meson) can be represented as follows: 

da D _ da D IF + da r ^ c 
dx dx dx ' 

da MF 

where -£ is the differential cross section production for D-meson due to c- 
quark fragmentation (that is described by the equation ([□]) with the modified 
fragmentation function of D MF (z, M c5 )), while is the differential cross 
section production for D-meson resulted from c-quark recombination with 
valence quarks (see (|9|)). 

Like in the previous article [0], we assume that a mesonic state (cq) goes 
over into a vector My or a pseudoscalar Mps meson with the probability 
proportional to the spin factor: 

M PS : My = 1 : 3. (20) 

3 COMPARISON OF THE MODEL PREDIC- 
TIONS WITH THE EXPERIMENTAL RE- 
SULTS IN 7T- BEAMS 

As it mentioned in the Introduction our article deals with the description 
of two types of the inclusive x-distributions, namely, the differential cross 
section for D meson production (i.e. 4|) and the asymmetry A{x). 

Fig. 3 presents the description of the differential distribution ^ for the 
reaction 

tt- N -> (D + D) X. 



(19) 
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The experimental data is summed with respect to all types of D meson. The 
beam energy equals E n = 250 GeV. 

As is seen from the figure our model (the modified fragmentation + re- 
combination) reproduces satisfactorily the experimental data. Note, also, 
that although the recombination contribution into the total cross section is 
rather small (< 10%), it plays a substantial role at high x (see Fig. 3). In 
addition, one may deduce from this figure that pure modified fragmentation 
can not provide the proper description of the inclusive spectrum in the whole 
kinematic region. 

Fig. 4 presents the description of the corresponding asymmetry (the 
leading particle effect). In this case, too, the considered model provides the 
description of the experimental data [[LT], [TTj . 

It should be particularly emphasized once more that the "hard" frag- 
mentation (i.e. . ftp = p c ) makes it possible the description of the differ- 
ential spectrum in the whole kinematic region but can not reproduce the 
x-dependence of the asymmetry. 

In the Table 1 we compare our model predictions with the experimental 
data on the total yields of the charmed mesons as well as with the predic- 
tions of the Lund model |12[] (the experimental results and the Lund-model 
predictions are taken from [|i~3|j ). 

One may deduce from the Table 1 that our evaluations for the ratios of the 
cross sections production of D mesons are in agreement with the experimental 
results within the experimental errors. Remind, that we obtain these results 
by consideration of the two types of the hadronization of charmed c quarks. 
The behavior of these processes is determined by two parameters R and M , 
as well as the distribution functions of the partons in the initial hadrons. 

The decisive test of the considered model should be the comparison of the 
theoretical predictions with the experimental results in and S~ beams. 
The valence quarks distributions in these hadrons differ substantially from 
those in 7r ± and p beams (see below). Consequently, one should expect 
the different contributions due to recombination mechanism in the inclusive 
spectra of charmed hadrons. 
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4 CHARM PRODUCTION IN THE BEAMS 
OF CHARGED K± MESONS 

From the view of parton model the difference of K ± mesons from mesons 
not only resides in the replacement of the valence d-quark by the strange 
valence s-quark. The distribution function of valence quarks in mesons 
should change substantially, as well. 

The simplest (without scaling violation) parameterization of the distri- 
bution functions of the valence q\ quark in the meson of M(qiq 2 ) with the 
valence quarks of q% and q 2 has the form as follows ||: 

v M( qi q 2 )( T \ _ r ( 2 + 7o ~ «i ~ «2) -ax,-, _ T yro-a 2 / 21 \ 

91 { ) " r(i - ai )r(i + 70 - a 2 ) 1 j ' { } 

where aij is the intercept of the leading Regge-trajectory for the quark, 
while 70 is certain parameter. 

The coefficient in the above equation is determined by the normalization 
condition: 



o 



V q M (x)dx = I. 



The choice of the parameter of 70 in ([H]) is governed by asymptotic 
behavior of the structure functions at x — > 1. The quark counting rules 
predicts the value as follows: 

70 - a 2 = I. 
From the well-known asymptotic for 7r-meson, 

Jx 



if follows that 70 = §• Taking into account that a u = ad = 1/2, while a s ~ 
(see (|8|)), we obtain the following distributions of the valence u and s quarks 
in K ± meson: 

V U K ~ 4=(l-*) 3/2 , (22) 
V x 

V S K ~ (1-x) 1 . (23) 
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From the form of the distributions of ( p2| ) and ( p3[ ) it is evident that the 
valence s-quark in the K meson is much "harder" than u quark: 



< xf > = 0.33, 

< xf > = 0.166. 

From the given parameterization of the K meson structure functions it follows 
that the total momentum fraction carried away by the valence quarks is equal 
to: 

< xf >=< xf v > + < < >= 0.5. (24) 
Such a magnitude should be compared to the similar value for the 7r meson: 

< < >=< xl> + < xl v >= 0.4. (25) 

We assume further that the distribution of the gluons in tt^ and K mesons 
are identical in form. The previous analysis |14[ of the evolution of structure 



functions of 7r ± and K mesons, where the evolution starts from the different 
distributions of valence quarks, provides an argument in favor of such an 
assumption. Other arguments in support of such an assumption are the 
identical form of the spectra of charmed mesons produced in 7r ± and 
beams [|1(] . 



As mentioned above, the form of the distributions of the rest sea partons 
in i^-meson coincides with the form of corresponding parton distributions 
in tt^ mesons: 

Here, the parameter e takes into account the variation of the momentum 
fraction carried away by the valence quarks in K meson with respect to to 
analogous value in n ± meson: 



1- < xf > 
1- < xl > 



0.8. 

Due to the different distributions of fl2"2] ) and (^) we should expect the 
different x-distributions of D(cu) and D s (cs) mesons, produced in the K ± 
beams (that resulted from the recombination with valence quarks). Two- 
particle distributions of the partons in K meson, needed for such a calcula- 
tion, can be easily evaluated with the help of equation (pT|) . Their explicit 
form is presented in the Appendix 2. 
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Fig. 5 presents our model predictions for the leading particle effect in K~- 
meson beam at the energy of E K = 250 GeV. As expected, our model predicts 
the different x-dependence of the asymmetry for D and D s mesons. For the 
strange-charmed _D s -mesons the leading particle effect is more pronounced as 
compared with usual charmed mesons. The observation of such a difference 
would be evidence in favor of the considered model. 

Unfortunately, only the integral asymmetry magnitude is measured in the 
experiment |15| : 



Ak P (D s ) = 0.25 ±0.11, (26) 

that should be compared with our prediction: 

A% eor {D s ) = 0.29. (27) 

As is seen our theoretical estimate is in agreement with the experimental 
value. 



5 CONCLUSION 

This article presents the "improved" model of charmed quark hadronization. 
Our model provides the consistent description of the inclusive differential 
spectra of D mesons, produced in n~p collisions. Further step forward in 
the understanding of a hadronization mechanism of heavy quarks would be 
related with a consideration of the processes of charmed particle production 
in K and E interactions. It caused, in particular, the different distributions 
of valence quarks in K and £ hadrons as compared with tt and p beams. The 
next step of such an investigation should be a consideration of the process 
of charmed baryons production, where the diquarks from the initial hadrons 
play a substantial role. 
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APPENDIX 1. 

In order to deduce the parameterization of a(M c5 ), which takes into ac- 
count the conditions (|TS|), i.e. a{2m c ) = oo and a(M ) = a c , we consider 
the expression for the first moment of \i from the fragmentation function of 
D MF (z,M c5 ) from fllT 



^(M cE ) = / zD(z,M c5 )dz 
Jo 



3-a(M c5 )' 
The desired expression for a(M cg ) equals: 

,., , 1 - 3<i(M a ) 
Q(M "- ) = 1 - MAfc) ' 

We assume the following (Q CD-motivated) dependence for /i(M cg ) 



In g 



where d ~ 0.464 is the parameter, similar to anomalous dimension. 
The new parameter q is expressed though M as follows: 



where 



M \ — 



2m, 



i 1 — Ot 

i/ = /i(M )~3 and /i(M ) 



3 — a r 
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APPENDIX 2. 

For completeness sake we present here the explicit form for two-particle 
distribution functions fy^xy^Xi) in ti ± and K mesons. 

Note, that two-particle parton distributions (like the single-particle func- 
tions) can not be theoretically evaluated. Therefore, we use the simplest 
phenomenological expression that takes into account the momentum conser- 
vation, the (1 — Xi) n dependence of sea partons, as well as the normalization 
conditions as follows (see for the details): 

(l-si) 

f Vi {x v ,xi)dx v = fifa), 

where fi(xi) is the single-particle distribution of an i-type parton in the 
hadron h. 

For two valence quarks such a distribution has the form ||: 

/Wfzi xo) - r ( 2 + 7o-«i-«2) Q1 a2( , 70 _i 
fvv[XuX2) ~T(l- ai )T(l-a 2 )T( l0 ) Xl X2 [L Xl X2} • 

For the case of one valence and one sea parton the corresponding distribution 
equals: 

fvj{x v , Xj) = Nj—— yFT-i i II \ x v Vx j (1 ~~ x v ~ x j) "(1 — Xj) , 
1 (1 — a v )L (1 + n v ) 

where Nj is the corresponding normalization factor of a sea parton, n v = 
7o — a i — 0C2 + (x v , and k = rij — 1 — 70 + aci + a.%. 

The parameters of the function fvj{%i->%2) for tt and K mesons are pre- 
sented in the Table 2. 
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Table 1. The ratios of the total yields (at x > 0) of the charmed mesons 
in the considered model, the experimental results [O] and the Lund-model 



predictions (taken also from |T3||). 





D++D- 




D~ 


D° 




D°+D° 


D°+D°+D++D- 


D+ 


w 


our model 


0.332 


0.102 


1.16 


1.0 


experiment 


0.416 ±0.016 


0.129 ±0.012 


1.35 ±0.05 


0.93 ±0.03 


Lund 


0.472 


0.077 


2.25 


1.09 



Table 2. Values of the parameters appearing in the functions fvj{xi, x 2 ) = 
Ax^ ai x 2 a2 (1 - x x - x 2 ) n (l - x 2 ) k . 



ix^~ meson 


A 


«i 


a 2 


n v 


k 


partons 












%i v d v 


0.477 


0.5 


0.5 


0.5 





u v g 


1.50 


0.5 


1.0 


1.0 


1.5 


U v (u, d?) sea 


0.090 


0.5 


1.0 


1.0 


3.5 




0.045 


0.5 


1.0 


1.0 


3.5 


K meson 


A 


«i 


a 2 


n v 


k 


partons 












TlySy 


1.27 


0.5 


0.0 


0.5 







1.34 


0.5 


1.0 


1.0 


1.5 


s v g 


3.22 


0.0 


1.0 


1.5 


1.5 


U v [Uj d)sea 


0.08 


0.5 


1.0 


1.0 


3.5 


y>vSsea 


0.04 


0.5 


1.0 


1.0 


3.5 


S v (u, d)sea 


0.192 


0.0 


1.0 


1.5 


3.5 




0.096 


0.0 


1.0 


1.5 


3.5 
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d a / d x 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

X 



Figure 1: Differential distributions jj~ for charmed c quarks for all 
values of M c5 (the upper histogram) and for M c5 > Mo = 10 GeV (the 
lower histogram). The experimental points correspond to charmed 
particles yield, summed over all types of D and D mesons (the reaction 
of tt'N collisions at = 250 GeV @.) 
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Figure 2: The description of the asymmetry A(x) (the leading par- 
ticle effect) in ir~p collisions [10, 11] in the mechanism of the "hard" 



fragmentation (see (|l4D). The different histograms correspond to dif- 
ferent values of the parameter Rq for the recombination mechanism 
(see (§)). The value of Ro = 3.0 corresponds to the upper histogram 
and so on. 
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d a/d x 




Figure 3: Differential distributions ^ for the energy of E n = 
250 GeV. The experimental data are taken from jL"0|]. The dotted 
(dashed) histogram corresponds to the recombination (fragmentation) 
contribution. The solid histogram represents their sum. The cross sec- 
tions are presented in /ib. 
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Figure 4: The description of the asymmetry A(x) in ir~p col- 
lisions [[n], 11] in the mechanism of the modified fragmentation 
(see (|19)) with the help of fragmentation function from (|l~7|)). 
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Figure 5: The predictions of the modified fragmentation mechanism 
for the asymmetry A(x) in K~p collisions at Ek = 250 GeV. The 
histograms correspond to the ratios of: 1) D s {cs) to D s (cs) mesons, 
2) all charmed D + D s mesons, 3) D s (cs) to all charmed D + D s 
mesons. 
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